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Abstract 

No verbal explanation can indicate a direction in space or the orien- 
tation of a coordinate system. Only material objects can do it. In this 
article we consider the use of a set of spin-| particles in an entangled 
state for indicating a direction, or a hydrogen atom in a Rydberg state 
for transmitting a Cartesian frame. Optimal strategies are derived for 
the emission and detection of the quantum signals. 



1 Formulation of the problem 

The role of a dictionary is to define unknown words by means of known ones. 
However there are terms, like left or right, which cannot be explained in this 
way. In the absence of a formal definition, material objects must be used to 
illustrate these terms: for example, we may say that the human liver is on the 
right side. Likewise, the sign of helicity may be referred to the DNA structure, 
or to the properties of weak interactions. 

Here we consider some cases where information cannot be explained ver- 
bally. The simplest one is when the emitter (conventionally called Alice) wants 
to indicate to the receiver (Bob) a direction in space. If they have a common 
coordinate system to which they can refer, or if they can create one by observ- 
ing distant fixed stars, Alice simply communicates to Bob the components of a 
unit vector n along that direction, or its spherical coordinates 9 and 0. But if 
no common coordinate system has been established, all she can do is to send a 
real physical object, such as a gyroscope, whose orientation is deemed stable. 

2 Quantum transmission of a spatial direction 

In the quantum world, the role of the gyroscope is played by a system with large 
spin. For example, Alice can send angular momentum eigenstates satisfying 
n • J|-0) = This is essentially the solution proposed by Massar and 

PopescuEl who took N parallel spins, polarized along n. The fidelity of the 
transmission is usually defined as 

F=(cos2(x/2)) = (l + (cosx))/2, (1) 
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where x is the angle between the true n and the direction indicated by Bob's 
measurement. The physical meaning of F is that \ — F = (sin^(x/2)) is t 
mean square error of the measurement, if the error is defined as sin(x/2) 
The experimenter's aim, minimizing the mean square error, is the same as 
maximizing fidelity. We can of course define "error" in a different way, and then 
fidelity becomes a different function of x Emd optimization leads to different 
results H. Here, we shall take Eq. (1) as the definition of fidelity. 

Massar and Popescu showed that for parallel spins, \ — F = 1/{N + 2). It 
then came as a surprise that for N = 2, parallel spins were not the optiraal sig- 
nal, and a slightly higher fidelity resulted from the use of opposite spins □. The 
intuitive reason given for this result was the use of a larger Hilbert space (four 
dimensions instead of three). This raises the question what is the most efficient 
signal state for N spins, whose Hilbert space has 2^ dimensions. Will F ap- 
proach 1 exponentially? Actually, the optimal result is a quadratic approach, 
as shown below. 

Our first task is to devise Bob's measuring method, whc^e.mathematical 
representation is a positive operator-valued measure (POVM)B'EI. For any unit 
vector n, not necessarily Alice's direction, let |j, TO(n)) = | j, m(0, 0)) denote 
the coherent angular momentum stateQ that satisfies 



and 



We then haveL 



J^\j,m{n))^j{j + l)\j,m{n)), (2) 
n- J\j,m{n)) ^m\j,m{n)). (3) 



(2j + l)Jde^ I j, m{e, 0)) (j, m{e, cj>) 1 = 1^, (4) 

where 

dg^ := sin 9 d9 d4'/ Ait, (5) 

and Ij is the projection operator over the (2j-|-l)-dimensional subspace spanned 
by the vectors |j, m{9, (j))). If TV = 2, so that j is or 1, the two resulting sub- 
spaces span the whole 4-dimensional Hilbert space. For higher N, all the 
rotation group representations with j < N/2 occur more than once. We then 
have, if we take each j only once, from or i to A^/2, 



(N + ^r iN + l){N + 3) ^ 
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for even or odd N, respectively. For large N, the dimensionality of the acces- 
sible Hilbert space tends to and this appears to be the reason that the 
optimal result for 1 — _F is quadratic in iV, not exponential. 

We now turn to the construction of Bob's POVMH. Let p denote the 
initial state of the physical system that is measured. All these input states 
span a subspace of Hilbert space. Let 1 denote the projection operator on 
that subspace. A POVM is a set of positive operators i?^ which sum up to 
1. The index /x is just a label for the outcome of the measuring process. The 
probability of outcome n is tr(p£'^). In the present case, /i stands for the 
pair of angles 6(p that are indicated by Bob's measurement. If we want a high 
accuracy, these jputput angles should have many different values, spread over 
the unit sphere For example, the components of a continuous POVM, as in 
Eq. (^, are given by 

Ee^ = {2j + l)dg^ I j, m{e, 0)) (j, m((?, ^) |. (7) 

Such a POVM with m = j corresponds to the method of Massar and Pxipescu0. 
The choice m = j is not optimal. As shown by Gisin and PopescuQ for the 
case N — 2, signal states with opposite spins give a higher fidelity. With our 
present notations, these states are (|0, 0) + \l,0{n)))/\/2. They involve two 
values of j, but a single value of m, namely 0. 

One possibility to include several values of j in a POVM is to take a sum of 
expressions like (^). This brings no advantage, because a convex cppibination 
of POVMs cannot yield more information than the best one of thcmB. Optimal 
POVM components are usually assumed to have rank one. In the present case, 
each one of them should include all relevant j: 

Ee4,:=de^\e,(t>){e,4>l (8) 

where 

N/2 

1^,0) ^ y27TT|j,m(0, </))). (9) 

To verify that this is indeed a POVM, we note that in J Eg^ there are 
diagonal terms (2j + m{d, (t))){j, m{9, 0)|, which give Ij, owing to Eq. (^. 
The off-diagonal terms with ji ^ j2 vanish, as can be seen by taking their 
matrix elements between {ji,rni\ and 1^2, ^712) in the standard basis where Jz 
is diagonal. We haveli^ 

O2,m(0,0)|j2,m2) =I?HlW), (10) 
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with a similar (complex conjugate) expression for (ji, toiIji, to(0, 0)). The 
rotation matrices T) are explicitly given by 



(11) 



where the Euler angle i}} is related to an arbitrary phase which is implicit 
in the definition of |j, m(6', 0)). Note that a single value of m occurs in all 
the components of the vectors \9,(j)) in Eq. (||), so that the undefined ^phases 
giimi/i niutually cancel. It then follows from Eq. (4.6.1) of Edmonds t2l that 
all the off-diagonal matrix elements of jEg^ vanish, so that we indeed have a 

povmb. 

While Bob's optimal POVM is essentially unique in the Hilbert space that 
we have chosen, Alice's signal state, which is 



N/2 

contains unknown coefficients c,-. The latter are normalized, 



(12) 



N/2 

j=m 



(13) 



but still have to be optimized. 

The probability of detection of the pair of angles 6'0, indicated by the 
POVM component Eg^, is 



{A\Ee^\A) = de^ 



N/2 



CjV2j + 1 {j,m{e,(j))\j,m{n)) 



(14) 



We havet 



{j,mie,cf>)\j,min))=e^^S±ix), (15) 

where x is the angle between the directions n and Ocf), and the phase e*'' is 
related to the arbitrary phases which are implicit in the definitions of the state 
vectors in (15). The important point is that e^^ does not depend on j and 
therefore is eliminated when we take the absolute value of the sum in Eq. (p^. 
Explicitly, we have 

d(^L(x) = cos2'"(x/2) P^^yico^x), (16) 
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where P^'^^ [x) is a Jacobi polynomialQM. We shall write x — cos x for brevity, 
so that the fidelity is 

F={l + {x))/2. (17) 

Our problem is to find the coefficients Cj that maximize {x). Owing to ro- 
tational symmetry, we can assume that Alice's direction n points toward the 
z-axis, so that de^ can be replaced by dx/2 after having performed the inte- 
gration over (f). We thus obtain 



(x) = i / xdx 



-1 



N/2 
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(18) 



This integral can be evaluated explicitly by using the orthogonality and 
recurrence relations for Jacobi polynomials E3. The result is 

(x) = ^c*Cfe^j-fc, (19) 
where Ajk is a real symmetric matrix, whose only nonvanishing elements are 



^3 3 

and 



A„=m''/[j{j + l)], (20) 



= = if - m2)/jV4j2-l. (21) 

The optimal coefficients cj are the components of the eigenvector of Ajk that 
corresponds to the largest eigenvalue, and the latter is (x) itself. For m = 
(which is the optimal choice) and large N, we find that 

1 -F^ [2.4048/(iV-h3)]^ (22) 

where the numerator is the first zero of the Bessel function JqEI Tie right 
hand side ought to be compared to the result of Massar and Popescutl, which 
was 1/{N + 2). A detailed comparison oOjoth methods is illustrated by the 
figure which appears in a previous articlellil. 

3 Transmission of a Cartesiam frame 

Next, we consider the quantum transmission of a complete Cartesian frame. If 
many spins are available, a simple possibility would be for Alice to use half of 
them for indicating her x axis and the other half for her y axis. However, the 
two directions found by Bob may not then be exactly perpendicular, because 



5 



separate transmissions have independent errors due to Umited angular resolu- 
tion. Some adjustment will be needed to obtain Bob's best estimates for the 
X and y axes, before he can infer from them his guess of Alice's z direction. 
This method is not optimal, and it is obviously not possible to proceed in this 
way if a single quantum messenger is available. Here we shall show how a sin- 
gle hydrogen atom (formally, a spinless particle in a Coulomb potential) can 
transmit a complete frame. 

Consider the n-th energy level of that atom (a Rydberg state). Its de- 
generacy IS d = in? because the total angular momentum may take values 
j = 0, • • • , n — 1, and for each one of them m = —j, ■ ■ ■ ,j- Alice indicates her 
xyz axes by sending the atom in a state 

n-l j 

l^> =E E <'jrn\j,m), (23) 

j=0 m=-j 

with normalized coefhcients a~ that will be specified below. Bob then per- 
forms a covariant measurementllj in order to evaluate the Euler angles ip6(p that 
would rotate his own xyz axes into a position paralleLto Alice's axes. Bob's 
detectors (ideally, there is an infinite number of themE^I) have labels ipOtj) and 
the mathematical representation of his apparatus is as usual a POVM. In the 
present case, we have a resolution of identity by a set of positive operators: 

J d^g^ E{ipe(j)) = 1, (24) 

where d^g^ = sm0dtpd0d(f>/8'!T^ is the 5*0(3) Haar measure for Euler ang les0, 
and E{il}9(j)) — \ij}94)) {il:9(j)\. The vectors \'4'9(j)) will be specified below. The 
probability that the detector labelled is excited is given by 

P{'4^9dp) = {A\d^g^E{^p9(j,)\A) = d.^g^ \{A\-^9(j))\'' . (25) 

Our task is to construct vectors \4'd<p) such that Eq. (^) is satisfied (that is, 
the probabilities sum up to one) and Bob's expected error is minimal. 

Generalizing the method of the preceding section, we define a fiducial 
vector for Bob, 

n-l j 

\B) - v/27TT fe,,„ |j, m), (26) 

j=0 'in=-j 

where the coefficients bjm are normalized for each j separately: 
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E l^^™l' = l ^J- (27) 

For the preceding problem, a single value of m was used; here we need all the 
values. Note that Eq. ( p3| ) was written with Alice's notations (m is the angular 
momentum along her z axis), while Eq. ( p6| ) is written with Bob's notations 
(m refers to his z axis). This issue will be dealt with later. 
We now define 

\m) = u{m) \B), (28) 

where U {ipO(j)) is the unitary operator for a rotation by Euler angles ipdcf). Note 
that since \B) is a direct sum of vectors, one for each value of j, then likewise 
U{il)9(j}) is a direct sum with one term for each irreducible representation, 

j 

where the D^-''^ (tpOcf)) are the usual irreducible unitary rotation matrices0. To 
prove that Eq. ( p^ is satified, we note that its left hand side is invariant if 
multiplied by U{iJivp) on the left and U{iJivp)^ on the right, for any arbitrary 
Euler angles p,vp (because these unitary matrices represent group elements and 
therefore have the group multiplication properties) 13. It then follows from a 
generalization of Schur's lemma EZl that the left hand side of ( p4[ ) is a direct 
sum of unit matrices, owing to the presence of the factor (2j + 1) which is 
the dimensionality of the corresponding irreducible representation. Therefore 
Eq. ( p^ ) is satisfied. 

The detection probability in Eq. ( [25| ) can thus be written as P{'ip9(j)) = 
djjierj, \ {A\U[ip9(l))\B)\'^ . To compute this expression explicitly, we must use con- 
sistent notations for \A) and \B) — recall that Eq. ( ^3| ) was written in Alice's 
basis, and Eq. (^ ) in Bob's basis. It is easier to rewrite Alice's vector |^) in 
Bob's language. For this we have to introduce the Euler angles £,rjC, that rotate 
Bob's xyz axes into Alice's axes (that is, S^rjC, are the true, but unknown values 
of the angles ^9(j) sought by Bob) . The unitary matrix U {^ri() represents an 
active transformation of Bob's vectors into Alice's. Therefore, U{£,ri()'' is the 
passive transformation El from Bob's notations to those of Alice, and C/(^?7C) 
is the corresponding transformation from Alice's notations to Bob's. Written 
in Bob's notations, Alice's vector \A) becomes U{£^ri()\A) so that, in Eq. (p5|), 
{A\ becomes {A\U{^ri(y . Let us therefore define 

UiaP-f) = Ui^riC)^ U{tp9(j)). (30) 
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The Euler angles a(3^ have the effect of rotating Bob's Cartesian frame into 
his estimate of Ahce's frame, and then rotating back the result by the true 
rotation from Alice's to Bob's frame. That is, the angles a/37 indicate Bob's 
measurement error. The probability of that error is 



P(a/37) = \{A\U{ap-i)\B)\\ (31) 

where dap-y = sin (Hdadpd'y /Stt'^ . Note that in the above equation \A) is written 
with Alice's notations as in (p3|), and \B) with Bob's notations as in (p^). 

Of course Bob cannot know the values of af3j. His measurement only 
yields some value for ip9(j). The following calculation that employs afHj has 
the sole purpose of estimating the expected accuracy of the transmission (which 
does not depend on the result tjjOcj) owing to rotational symmetry). 

We must now choose a suitable quantitative criterion for that accuracy. 
When a single direction is considered, it is convenient to define the errorH as 
sin(a'/2), where lu is the angle between the true direction and the one estimated 
by Bob. The mean square error is 



(sin2(w/2)> = (1 - (cosw))/2 = 1 - F, (32) 

where F is the fidelity, as defined above. When we consider a Cartesian frame, 
we likewise define fidelities for each axis. Note that coswfe (for the fc-th axis) 
is given by the corresponding diagcmal element of the orthogonal (classical) 
rotation matrix. Explicitly, we haveliS 



coswz = cos/3, (33) 

and 

coswa; + cosojy = (1 + cos/?) cos(q; + 7), (34) 
whence, by Euler's theorem. 



cos ujx + cos ujy + cos [^2 = 1 + 2 cos n, (35) 

where has a simple physical meaning: it is the angle for carrying one frame 
into the other by a single rotation. 

The expectation values of the above expressions are obtained with the help 
of Eq. dH): 

(/M7)) = J do.p, |(^|C/(a/37)|i3)P /(a/37), (36) 
where, explicitly. 
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{A\U{a(3^)\B) = a*^b,r{j,m\V^^\aP^)\j,r). (37) 

The unitary irreducible rotation matrices T>^^\af3^) have components0 

(j, m|I?W(a/?7)|j,r) =e'('""+'-^)dW(/3), (38) 

where the dWriP) can be expressed m terms of Jacobi polynomials. Collecting 
all these terms, we finally obtain, after many tedious analytical integrations 
over products of Jacobi polynomials £3, 

(/(a/37)) = X! fjkmnrs bjr akn bis, (39) 

where the numerical coefficients fjkmnrs depend on our choice of f{af3^) in 
Eqs. (^3^351). The problem is to optimize the components aj,„ (normalized 
to 1), and bjm satisfying the constraints (p7|), so as to maximize the above 
expression. For further use, it is convenient to define a matrix 

^ jm^kn — ^ ^ fjkmnrs bjr bf^^i (^^) 



r,s 



SO that 

(/(a/37)) = E ^om.kn a*„ afc„ = {A\M\A). (41) 

First consider a simple case: to transfer only the z axis, we wish to maxi- 
mize (cos /3) . An explicit calculation yields 

fjkmnrs — ^mn ^rs 9jk- (^^) 

The matrix gjk which is defined by the above equation has nonvanishing ele- 
ments 

5,j=W[jXi + l)], (43) 

and 

_.,a^SEZl^ (44) 



The 5mn term in (42) implies that, for any choice of Bob's fiducial vector \B), 
the matrix M in ([4l|) is block-diagonal, with one block for each value of m. The 
optimization of Alice's signal results from the highest eigenvalue of that matrix. 
This is the highest eigenvalue of one of the blocks, so that a single value of m 
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is actually needed. A similar (slightly more complicated) argument applies if 
Alice's vector is given and we optimize Bob's fiducial vector. This result proves 
the correctness of the intuitive assumption that was made in our solution of the 
first problem, for which a single value of m was used. It was then found that 
when m = (the optimal value) the expected error asymptotically behaves as 
1.446/d, where d is the effective number of Hilbert space dimensions. In the 
present case, d = (jmax + 1)^, which is the degeneracy of the n-th energy level. 

If we want to transfer two axes, we use Eq. (Q) and calculate the matrix 
elements for ((1 + cos (3) cos (a + 7)). (It is curious that they are simpler than 
those for (cos (a + 7)) alone.) We obtain 

f jkrnnrs — ^m,n— 1 ^r,s— 1 hj}^ ^ ^s,r— 1 h}^j^ (^^) 

where the nonvanishing elements of hjk are 



[(j -n+ +n)(j-s^ l)(j + .s)]V2 ^^^^ 

(47) 
(48) 

Note that the hj^ matrix, whose elements depend on n and s, is not symmetric 
(while was). This is because it comes from the operator e*^"^'''-' which is 



"■3 J ~ 


2i(j 4 


1) 








- s + 


1)(J-.)]V2 




2j{Af - 


1)1/2 




[(j + n - 




+ s - 


l)(j+.s)]l/2 




2j{Af - 


1)1/2 





not Hermitian. However the two terms of (45) together, which corresponds to 
cos(a + 7), have all the symmetries required by the other terms a*„j hjr dkn b'^s 
in Eq. (|39|). Finally, if we wish to optimize directly the three Cartesian axes 
(without losing accuracy by inferring z from the approximate knowledge of x 
and y) we use all the terms of (^5|), that is, both those of ( p2| ) and of (^5|). 

It now remains to find the vectors \A) and \B) that mipipiize the transmis- 
sion error. For small values of j, we used Powell's methodlia without imposing 
any restrictions on \A) and \B) other than their normalization conditions. As 
intuitively expected, we found that the optimal vectors satisfy 

/ N-l/2 
bj,n = aj,n ^ IfljnP , Vj. (49) 

This means that Bob's vector should look as much as possible as Alice's signal, 
subject to the restrictions imposed by the constraint (E|). 
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Taking this property for granted is the key to a more efficient optimization 
method, as follows: assume any bjm, so that the bilinear form ( pO| ) is known. 
Find its highest eigenvalue and the corresponding eigenvector ajm- From the 
latter, get new components bjm by means of (49), and repeat the process until 
it converges (actually, a few iterations are enough). Quantitative results are 
displayed in the figure of a preceding articles. That figure shows that there is 
little advantage in optimizing only two axes, if for any reason the third axis is 
deemed leas important. If the three axes are simultaneously optimized, it can 
be shownEj that the mean square error tends asymptotically to l/Vd. 

It is not surprising that this result is weaker than the one for a single 
axis, which was 1.446/d. The obvious reason is that we are now transmitting 
a three-dimensional rotation operation that can be applied to any number 
of directions, not only to three orthogonal axes. Indeed, consider any set 
of unit vectors e^, where m = 1,2,3, and is a label for identifying the 
vectors. Let be a positive weight factor attached to each vector, indicating 
its importance. Let R{a(3j) be the classical orthogonal rotation matrix ta for 
Euler angles a(3j. Then the cosine of the angle between Bob's estimate of 
and the true direction of that vector is 

cosw^ = ^i?™„(a/37)e^„e^. (50) 
With the same notations as before, we have 

(/(a/37)) = ^^tJ. (COSCJ^) = ^{Rmn{a(3^)) Cmn, (51) 

where 

Cmn=Y.'^t^<n<- (52) 

This is a positive matrix which depends only on the geometry of the set of 
vectors whose transmission is requested. We can now diagonalize Cmn and 
write it in terms of three orthogonal vectors, possibly with different weights. 
Therefore, no essentially new features follow from considering more than three 
directions. 

Finally, we note that all the above calculations, as well as those in preced- 
ing works, assume that Alice and Bob have coordinate frames with the same 
chirality (this can be checked locally by using weak interactions). If the chi- 
ralities are opposite, then all the directions inferred by Bob should be reversed 
(because directions are polar vectors while spins are axial vectors). 
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In summary, we have shown that a single structureless quantum system (a 
point mass in a Coulomb potential) can transfer information on the orienta- 
tion of a Cartesian coordinate system with arbitrary accuracy. At first sight, 
this conclusion seems surprising. No spherically symmetric classical object can 
achieve this result. Only those having an asymmetric internal structure, such 
as an asymmetric rigid body, can reliably transmit a Cartesian frame. However, 
a classical Kepler orbit has two vectorial constants of the motion: the angu- 
lar momentum and the Lenz vector 113. The Hamiltonian itself is spherically 
symmetric, but each elliptic orbit (each solution of the equations of motion) 
defines a unique Cartesian coordinate system. From the correspondence prin- 
ciple which is valid for large quantum numbers, we expect that there are sets 
of coefficients aj„ such that the wave function is concentrated in the vicinity of 
a classical elliptic orbit, and thereby defines a Cartesian system. This problem 
is under current investigation. 

4 Covariant quantum measurements may not be optimal 

In the preceding section, we discussed the transmission of a complete Cartesian 
frame by a single quantum messenger, namely a hydrogen atom (formally, a 
spinless particle in a Coulomb piotential). A similar calculation was done by 
Bagan, Baig, and Muhoz-Tapia Ell (hereafter BBM), who were able to reach 
much higher values of j by sophisticated analytical methods. There is an 
essential difference between our work and that of BBM. We considered a single 
system, while BBM took N spins, and one irreducible representation for each 
value j of the total angular momentum. The maximum value is jmax = N/2, 
and then the mathematics are the same as for our Rydberg state, with jmax = 
n — 1. However, as explained above, if there are N spins that can be sent 
independently, there is a better method. Alice can use half of them to indicate 
one axis, and the other half for another axis. This method is far more accurate, 
especially if N is large. From Eq. (|2^), the mean square error for each axis 
is 5.783/(iV/2)2 = rather than A/SN which is the resuh with the 

method used by BBMlII. Similar results hold even for low values of N. Why 
is there such a discrepancy? 

In all the works that were mentioned above, and in many other ^imilar 
ones, it was assumed that Holevo's method of covariant measurements 113 gave 
optimal results. That method considers the case where Alice's signals are the 
orbit of a group Q, with elements g. Namely, if \A) is one of the signals, 
the others are \ Ag) — U{g)\A), where U{g) is a unitary representation of the 
group element g. The problem is to find optimal quantum states for Alice's 
signals and Bob's detectors. Originally, Holevo considered only irreducible 
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representations- -Now we know that in some cases reducible representations 
are preferablellj'Eil. However one then never needs to use more than one copy 
of each irreducible representation in the reducible one. 

We now turn our attention to Bob. The mathematical representation of 
his apparatus is as always a POVM, namely a resolution of identity by a set 
of positive operators: 

Y,Eh = \ (53) 

h 

where the label h indicates the outcome of Bob's experiment. In the case of 
covariant measurements, the labels h run over all the elements of the group Q 
(with a suitable adjustment of the notation in the case of continuous groups). 
Then the probability that Bob's apparatus indicates group element h when 
Alice sent a signal \Ag) is 

P{h\g) = {Ag\Ey,\Ag). (54) 

The method of covariant measurements further assumes that can be written 
as 

Eh^\Bh){Bhl (55) 

where 

\Bh) = U{h)\B). (56) 

Here, \B) is a fiducial vector for Bob (which has to be optimized) and U{h) is 
a representation (possibly a direct sum of irreducible representations) of the 
same group Q that Alice is using. 

All this seems quite reasonable (and this indeed usually works well) but, as 
the above example shows, this may not be the optimal method. In the above 
example, Alice's signals \Ag), for all possible positions of her axes, are SO{?>) 
rotations of a fiducial state \A) with j = 0, ...,n — 1. On the other hand Bob 
uses two separate POVMs, each one testing only one half of the spins. Each 
one of these POVMs also involves 5*0(3), but with lower values of the maxiraal 
j. Further discussion of this topic will appear in a forthcoming publication E3. 

Work by AP was supported by the Gerard Swope Fund and the Fund for 
Encouragement of Research. PFS was supported by a grant from the Technion 
Graduate School. 
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